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^3 ' We study by means of numerical simulations the velocity reversal model, a one-dimensional mechanical 

Cn ' model of heat transport introduced in 1985 by laniro and Lebowitz. Our numerical results indicate that 

this model, although it does not conserve momentum, exhibits an anomalous Fourier's law similar to the 
one s previously observed in momentum-conserving models. This is contrary to what is obtained from the 
solution of the Boltzmann equation (BE) for this system. The pair correlation velocity field also looks very 

O^ , different from the correlations usually seen in diffusive systems, and shares some similarity with those of 

OA ■ momentum-conserving heat transport models. 
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Introduction 



Understanding the steady state of a system in contact with two heat baths at unequal temperatures is a central question to 

j^ I the theory of non-equilibrium systems|l!]. With very few exceptions (such as zero range processes^ l3|l), non equilibrium 

Cj ■ systems exhibit long range correlations in their steady state ||3-01. These correlations have been studied theoretically by various 

I , approaches and measured experimentally. 

'^ ' For a one dimensional diffusive system of size L, for instance a stochastic lattice gas, these long range correlations take the 

G| following scaling formlMllsl-UOll: 
O . 

^: {A{n)Air,)...A{r„)), = -^F„ (^, ...^) , (1) 

T-H ■ 

^ • (where A{r) is an observable at position r such as the density or the energy) when the distances between the positions r^ are 

00 ' macroscopic (i.e. L ^ \ri — rj\ ^ 1). 

0^ . The macroscopic fluctuation theory developed by Bertini et al. lllll[l2ll allows one to write down the general equations satisfied 

C ' by the scaling functions i^„[|8|]. These partial differential equations are usually difficult to solve. The expressions of the Fn are 

^^ , however known in a number of examples iRl [l3l [l4ll . where they have been obtained either from exact solutions of microscopic 

t^^ ■ models or by integrating the partial differential equations derived from the macroscopic fluctuation theory. Diffusive systems 

^^ ', are also known to satisfy Fourier's law, meaning that, for large system sizes, the steady state flux (J) of energy through a system 

^— ^ ■ of size L scales like L^^: 



(JL)--^G(r,,Tb), (2) 



;h ' where G{Ta,Tij) is a function of the temperatures Ta and Tb of the heat baths at the two extremities of the system, which 



maintain it out of equilibrium; the function G{Ta, Tb) vanishes linearly with the difference Ta — Tb so that the flux becomes a 
gradient. 

Over the last 15 vears lll5l - l23ll it has been realized on the basis of numerical simulations that one-dimensional mechanical 
systems which conserve momentum do not satisfy Fourier's law. Instead, they exhibit a power law decay, called anomalous 
Fourier's law, of the average flux J with system size: 

{Jl) = L''-'G{Ta,Tb) . (3) 

So far, the exponent a has not been determined analytically for any microscopic model with non-quadratic interactions; however, 
numerical ll24i 12511 and analvtic 11191 I26l - t28ll calculations, based on mode coupling theory or on other approaches 11291 130(1 . indicate 
that, depending on the type of the non-linearities and on the accuracy of the simulations, a can take values ranging between 0.2 
and 0.5. The main two systems for which this anomalous heat conduction has been observed are; 



• 



• 



the Fermi-Pasta-Ulam model, a chain of N harmonic oscillators with an additional cubic (FPU-a) or quartic (FPU-/3) 
interaction potential llisl \3Ul [3111 . For the latter case, most estimates indicate that a '^ 0.25; 



a one-dimensional gas of N hard point particles with elastic collisions ll25ll (when the particles are identical, the collisions 
simply exchange the velocities of the incoming particles: hence the system has the same transport properties as an ideal 
gas, i.e. a ballistic transport with a = 1). For the collisions to be non-trivial, one choice is to consider particles of 
alternating masses 1 and 7712 7^ 1 (2-mass model). In this case, it is commonly found that a ~ 0.33 12411 . 

In comparison with diffusive systems, for which ^ holds, the main feature of these models is that their dynamics conserve 
momentum, and it is believed that this property plays a key role in the divergence of the heat conductivity |32:]. Note, however, 
that there is no proof for this behavior: the argument of [32] is not correct, although its conclusion is consistent with observations. 
In this article, we consider a simpler model originally introduced in [33i|, the velocity reversal model, in order to better 
understand which features could be at the origin of the anomaly. This model is, like the ID hard-particle gas, a system of N 
free-moving particles undergoing collisions. The collisions are, however, not elastic, but instead are given by the following 
simple rule: 

• when two particles collide with velocities of opposite signs, Vi+i < < Vi, their two velocities get reversed ( Ui ^' —Vi, 
Vi+i —J' —Vi+i at the collision); 



• 



when they collide with velocities of the same sign, the particles simply pass each other ( Wi — !> w^+i and w^+i — !> Vi if the 
particles are kept ordered from left to right). 



In contrast with the Fermi-Pasta-Ulam chain and the hard-particle gas, these dynamics do not conserve momentum. It should 
be noted, however, that the absolute values of the velocities themselves are conserved, leading in general to more conserved 
quantities than in standard momentum-conserving models. One should also note that, when the velocities are reversed at inde- 
pendent exponential times instead of at collisions, the system becomes diffusive ll33ll . 

Here, the N particles of the velocity reversal model are in a one-dimensional box of length L — N between two boundaries 
which play the role of heat baths. Whenever a particle hits a boundary, its velocity is changed as if the particle was reinjected 
instantaneously, with its velocity thermalized by the reservoir The reservoirs are thus described by the velocity p.d.f. of these 
"reinjected particles", Pa{v) and Pb{v). An obvious choice is to consider MaxwelUan reservoirs at temperatures Ta and T^, with 
velocity densities 

Pa(«)=e(f)^e-^andp,(7;)^0(-7;)Me-4,with0(7;) = |jJ°;;;>[J . (4) 

We will study below this "Maxwellian case" for T^ = 4 and T^ = 1. For a system which thermalizes well, one would 
expect other choices of these reservoir velocity densities to only affect small regions near the boundaries, but not to influence 
the macroscopic behavior of the system. For the velocity reversal model, however, the conservation of the absolute values of 
the particle velocities prevents thermalization from taking place: for instance, one can easily see that, for another choice of 
reservoirs, where particles are always reinjected with velocities w„ from the left reservoir and «(, for the right reservoir, 

Pa{v) = 5{v - Va) and pb{v) = 5{v + vt) , (5) 

the only possible velocities inside the system are ztva and ±Vb- We will study this "two-speed case" with Va = 2 and vt ^ 1. 

In a first part, we show, by studying the steady-state current of systems with 125 < iV < 8000 particles, that, both in the 
Maxwellian and the two-speed case, the velocity-reversal model exhibits anomalous Fourier's law, in contrast to the prediction 
of the associated Boltzmann equation ll33ll . 

We then present, in a second part, measurements of the steady-state correlation functions of the momentum density for 
Maxwellian and two-speed systems of 100 < N < 400 particles, as well as for a 2-mass hard-particle gas with mi = 1 and 
1712 = 1.6 for comparison: they indicate a size dependence very different from ([T]). 

I. ANOMALOUS FOURIER'S LAW FOR THE VELOCITY REVERSAL MODEL 

The main advantage of the velocity reversal model described above, compared to the hard-particle gas for instance, is that 
its associated Boltzmann equation can be solved analytically, thanks to the simpler form of its collision term. By solving this 
Boltzmann equation in the steady state with the appropriate boundary conditions, the average current between the reservoirs 
can be predicted, as done in 13311. The simplest case is the two-speed system ^ with velocities Va and v^, for which the 
Boltzmann equation relates the four densities f±a{x, t) and f±b{x, t) of the particles with respective velocities ±.Va and iff,, 
with Q < X < L — N the space coordinate : 



(dt + Vadx)f+a = {Va + Vb)if+bf-a - f-bf+a) = M^i i 
{dt - Vadx)f-a = -A{x,t) ; 
{dt + vtd^)f+t = -A{x, t) ; 
{dt -Vhd^)f^h ^ A{x,t). 
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Figure 1: Steady-state energy current Jm for the velocity-reversal model with 125 < N < 8000 particles, in the Maxwellian case with Ta — 4: 
and Tb = 1 (above) and in the two-speed case with Va ~ 2 and vt — 1 (below). In both cases, the current decreases as a non-integer power of 
the system size : over the range of sizes we considered, our data are well-fitted by Jn oc 7V~°'®^ in the Maxwellian case and Jjv oc jV""'®" 
in the two-speed case. 
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Figure 2: Average density (left) and energy (right) profiles for the two-speed velocity reversal model between reservoirs for Va 
and for system sizes TV — 2 ^ "^ 

/+a + /-a + /+6 + /-b and 



2 and v^ 



10'', 4 ■ 10^ and 8 ■ 10'' (continuous lines). According to the Boltzmann equation, these profiles (given by 

2 2 

'■fif+a + f-a) + ^{f+b + f-b) respectively) should be linear in x. 



Setting the time derivatives to 0, these equations can be solved ll33ll to determine the steady state of the two-speed system: 

J-a(Xj = ^J-zp^Tv+r'.-'- ~ ^) 
J+b{x) = „^+t,^ Iv+T^ 

yielding linear profiles for the four densities as well as a steady-state current J = ""vtv+T)" satisfying Fourier's law©. 
By a similar calculation, the Boltzmann equation also predicts that the energy current should also satisfy Fourier's law in the 
Maxwellian case. 

We studied numerically the Maxwellian dUi and the two-speed Q velocity reversal models for 125 < A^ < 8000 particles, 
with Ta = 4:, Ti, ~ 1 in the first case and Va = 2, Vb = 1 in the second. Because the true steady state is unknown, our 
measurements were performed by starting the system in an equilibrium configuration at T = 2 in the Maxwellian case, and by 
choosing the particle velocities uniformly at random among ±Va, ±w;, in the two-speed case. 

We then let the system evolve from that initial state, and sampled the time evolution of Jjy, piy{x) and En{x) : Jjv was 

measured by a time average of the instantaneous energy flux, Jn ~ ^ Si=i '^f' while pn{x) and En{x) were measured by 
counting the particles and energy within the boxes fc < .t < fc + 1 for k = 0, .., A^ — 1. We then took the long-time limits of 
these measured quantities to estimate their values in the steady state. 

As shown in figure [T] our data indicate that Jn decreases like a power law of A'^ for both choices of reservoirs for the system 
sizes we considered. They are consistent with (|3]l, with a ^ 0.3 for the Maxwellian reservoirs and a ^ 0.4 for the two-speed 
reservoirs. In the latter case, the density and energy profiles pN (x) and En {x) shown in figure|2]also differ noticeably from the 
linear profiles predicted by the Boltzmann approach, at least for the sizes we were able to study. 

II. LONG-RANGE CORRELATIONS IN ANOMALOUS SYSTEMS 

The main assumption in the derivation of the Boltzmann equation is that there are no correlations between the velocities of 
particles entering a collision in which their velocities are reversed. As it fails to predict the anomalous Fourier's law of figure[T] 
it is interesting to investigate the form of the steady-state correlations for the velocity reversal model. 

We measured numerically the steady-state two-point correlation functions of the momentum field, {p{x)p{y)) c, for the velocity 
reversal model in the Maxwellian (fig. O and in the two-speed (fig. |4|i cases, for systems of 100 to 400 particles. 

As in the previous section, we approximated the momentum density function pn{x) by its discretization over the boxes 
fc < a: < fc + 1 for < fc < A^. We took the stationary- state correlations {p{x)p{y)) c, which we measured for y ~ ^ and 
y = -2^, to be the long-time limits of their time evolutions for a system started in the arbitrary initial state described in the 
previous section. 

For diffusive systems, one would expect {p{x)p{y))c to follow a scaling of the form ([T]i: {p{x)p{y))c = jjP2{jf, -pf)' except 
for X = y with a variance {p{x)'^)c of order 1. For the velocity reversal models of figures|3]and|4l we observed {p{x)'^)c to be of 
order 1: however, when measuring {p{x)p{y)) c for y ~ N/A and y — iN/A, as in figures [3] and HI, we observed an additional 
anomalous gowth in {p{x)p{y)) c for x close to ?; in a region growing with system size, but becoming narrower on a macroscopic 
scale. 

The sign of {p{x)p{y))c seems to be always negative for y closer to the left (warmer) reservoir and positive for y closer to the 
right, colder reservoir. As shown in figures [3] and |4] where our data are multiplied by the system size A^, the correlations seem 
to decay slower than 1/A^ as A^ grows, over a region which seems to grow slower than the system size. 

In order to compare the behavior of these models with those of usual momentum-conserving models, we also measured the 
two-point correlation function of the momentum for a hard-particle gas of particles of alternating masses 1 and 7712 = 1.6 
between Maxwellian reservoirs at Tq = 4 and T;, = 1 : we observed a similar growth of the correlation functions at intermediate 
scales (fig. |5]). 

For this hard-particle gas (fig. |5]), due to a faster relaxation to the steady state, we were able to simulate much larger system 
sizes than for the velocity reversal model: at these sizes, 10'^ < A^ < 4 • lO'^, the anomalous growth of {pN{x)pN{y))c of figure 
|4]is compatible with a scaling of the form 

(pNix)pN{y))c OC -rT7r7rF2 



as shown in figure |6] The sizes we could reach for the velocity-reversal model of figures |3] and |4] did not exhibit such a clear 
scaling form. 

For the velocity-reversal model as well as for the 2-mass hard particle gas, we also measured the correlations of the energy 
field En{x), which exhibited a behavior very similar to the one shown above for the momentum field. 
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Figure 3: Two-point correlation function {p{x)p{y))c of the momentum field pn{x) for the velocity reversal model in the Maxwellian case 
^ with Ta = 4 and T^ = 1, for i/ = f (left) and y = 2^ (right). 
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Figure 4; Two-point correlation function {p{x)p{y))c of the momentum field pn{x) for the velocity reversal model in the two-speed case (|5j 
with Va = 2 and vt = 1, for y = ^ (left) and y = ^ (right). 



Conclusion 



The numerical simulations presented in this work show that the velocity-reversal model exhibits an anomalous Fourier's law, 
although it does not conserve momentum. Its density and energy profiles (fig. |2]) differ noticeably from those predicted by the 
Boltzmann equations. For the sizes we could achieve, however, they have not yet converged: this might question whether our 
data showing the anomalous Fourier's law (fig. [T]) correspond to the asymptotic regime. 

Within the system sizes we could reach, the two-point function of the momentum also looks very different from those of 
diffusive systems: they seem to be concentrated on a mesoscopic scale, much larger than the microscopic scale, but smaller than 
the system size. Within this mesoscopic region, the decay of the correlations seems to be anomalous as well, i.e. seems to be a 
non-integer power of the system size. 

We have checked that the pair correlations seem to have a rather similar behavior for the 2-mass hard particle gas. This 
leads us to believe that the shape of the correlation functions we have observed is probably another signature of the anomalous 
Fourier's law: it would therefore be interesting to see whether these anomalous correlations are also present for other models 
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Figure 5: Two-point function {p{x)p{y))c of the momentum field pn{x) for the 2-mass hard-particle gas with 7712 = 1.6 between reservoirs 



at Ta = 4 and U = 1, for y 



(left) and y - 



(right). 
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Figure 6: Rescaling of the two-point function of the momentum field {pn {x)pn {y)) c, for y — ^ (left) and y = ^ (right), for the 2-mass 
hard particle gas with m2 — 1.6 between reservoirs at temperatures Ta = 4 and Tt = 1. 

known to exhibit an anomalous heat conductivity, such as the Fermi-Pasta-Ulam chain. 

Looking at higher-order correlations would require major numerical efforts. They would however be very interesting to 
measure in order to guess what would replace ([T]) for systems exhibiting anomalous Fourier's law. 

The present work is purely numerical in nature. It would of course be very interesting to see what the existing theories 
explaining the anomalous Fourier's law would predict for the pair correlation functions we have measured. For diffusive systems, 
the scaling ([T]i of the correlations with system size is closely linked to the scaling form of the large deviation function of the 
system's density profileiS, lOJ: what such a scaling form would look like for systems exhibiting anomalous Fourier's law is also 
an interesting open issue. 
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